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Abstract — A Wiener filter can be interpreted as a cascade 
of a whitening- and an estimation filter. This paper gives a 
detailed investigates of the properties of these two filters. Then 
the practical consequences for the overall Wiener filter are 
ascertained. It is shown that if the given spectral densities are 
smooth (Hoelder continuous) functions, the resulting Wiener filter 
will always be stable and can be approximated arbitrarily well by 
a finite impulse response (FIR) filter. Moreover, the smoothness 
of the spectral densities characterizes how fast the FIR filter 
approximates the desired filter characteristic. H on the other 
hand the spectral densities are continuous but not smooth enough, 
the resulting Wiener filter may not be stable. 

I. Introduction 

The solution of the the optimal smoothing and prediction 
problem was published by Wiener, Hopf and Kolmogoroff [1], 
[2] in the 1940's. The well known formal solution of this 
problem (the causal Wiener filter) is determined from two 
given spectral densities with certain properties. The question 
is, how does these properties pass over to the Wiener filter, 
for instance: to the smoothness of the solution [3], as a 
mathematical attribute; to the stability of the filter, as a 
practical characteristic; or to the possibility to approximate 
the Wiener solution by an FIR filter, as an important practical 
and mathematical property. The answer of these questions is 
non-trivial because the operations, necessary to calculate the 
causal Wiener filter, are non-linear. 

In this paper the problem is completely solved for the most 
important smoothness class of spectral densities, namely for 
Hoelder continuous functions. Moreover, it is shown that the 
results hold not for continuous spectra, in general. In this 
paper, all investigations are done in the energy norm, but a 
similar treatment in the bounded-input bounded-output norm, 
which yields completely different results, was done in [4]. 

II. Notations & Problem Statement 
A. Some system theoretical aspects 

Consider a time invariant and time-discrete linear system S. 
If {hk} is its impulse response then the input-output relation 
y — Sx in the time domain is y„ — J^kL-oo ^kXn-k in which 
Xn is an input- and j/„ is the corresponding output sequence. 
If hk — for all fc < 0, the system is called causal. The 
so called S! -transform of an arbitrary sequence {xk} and the 



corresponding inverse i^-transform is defined by 

x{z) = XkZ^ and Xk ^ — i x(z)z~^''^'^^dz . (1) 

Note that the ^-transform is the usual ^"-transform where z 
is replaced by z^^. The ^-transform of an impulse response 
{hk} is called the transfer function h of the linear system and 
its frequency response /i(e~*'^) with lu G [— 7r,7r) is obtained 
by setting z — e~*" in its transfer function. From the definition 
of the f^-transform (QJ and from the Cauchy integral theorem, 
it is clear that {hk} is causal if the corresponding transfer 
function h{z) is an analytic function inside the unit disk V 
{z€C:\z\< 1}. 

In the i^-domain the input-output relation becomes 
y{e~"^) = h{e-'-'^)x{e-"^) for u £ [-7r,7r). The average 
energy of a sequence Xn is given by the square of its P- resp. 
L'^-norm ||x||2. The stability norm of a linear system S is 
defined by 

\\S\\e-= sup ||5a;||2 
Ibll2<i 

and it is well known that — \\h\\^ where \\h\\^ = 

sup^ |/i(e~''^)| is the usual supremum norm of the frequency 
response of S. For a causal h this norm can also be written 
ll'^lloo ~ supj^i^j^ Any linear system S is called 

(energy) stable if it has a finite energy norm < oo, 

i.e. if any input x with ||a;||2 < oo yields an output y with 

I|y|l2 < 

Now, all causal and stable transfer functions are collected 
in the Hardy space 7i°°(I?). This is the set of all complex 
functions which are analytic and bounded (with respect to 
II 11^) in v. Together with the point wise multiplication, 
TL°°{'D) is a Banach algebra. 

Consider now a transfer function hj^ of the form h]^{z) = 
'^k=o ^kZ^ ■ Practically, the corresponding Sn can be realized 
by an FIR filter of degree N . We say that a linear system S 
can be approximated uniformly by an FIR filter if there exists a 
Sn with transfer function hN such that the maximal difference 
between the output j/at = Snx of the FIR filter and the output 
y = Sx of the linear system becomes arbittarily small as N 
tends to infinity, i.e. if 

- ^Tvll^; = sup \\{h- hN)x\\^^ as ^ oo. 

Ikllt2<l 



It is known, that any h E 71°° {V) which is continuous in the 
closure of the unit disk T) can be approximated uniformly by 
an FIR filter. The sub-algebra of 7i°°{'D) which contains these 
functions is the disk algeba A{'D) C 7i°°(I?). This is the set 
of all functions analytical in T) and continuous in V. A{V) 
comprises all stable and causal transfer functions which can 
be approximated by FIR filters. 

Since we consider discrete-time systems, the frequency 
responses are defined on the unit circle, i.e. they are defined 
on [— TT, tt). Therefore, we denote by C[— 7r,7r) the set of 
all complex valued bounded and continuous functions on the 
interval [— 7r,7r) with /(— tt) = lim^^^Tr /(^) To control 
the smoothness of our functions, we use the so called Hoelder 
spaces (see e.g. [5]): Let be a domain in the complex plain. 
The set of all functions / for which 



l/L 1/(0)1 



|/(Z1)-/(Z2)| 

sup — . — < oo, 



is denoted by Ca(S^- This is the set of all / which satisfy a 
Hoelder condition of order a. In particular, we need the spaces 
Cq.[— TT, tt) and ^^(I?) with < a < 1. The latter is the subset 
of A(jy) which elements satisfy a Hoelder condition of order 
a. With the above norm they are Banach algebras. 

The set of all real valued functions /(w) defined on an 
interval with ^{lS) > for all w G O is denoted by C+(r2). 

B. Wiener filtering 

The Wiener filtering problem can be formulated as follows 
[6]: Given a wide sense stationary random process {?/„} with 
power spectrum $ as the input of a linear filter H. Let 
{xn} be the desired output of this filter and denote by the 
Fourier transform of the cross-correlation function between the 
sequences Xn and y„. We look for a causal H such that the 
output of this filter Xn = X^feLo Hky-n-k minimizes the mean 
square error (MSB) e = £[\xk — Xk^] to the desired sequence 
Xn, where £ [■] stands for expectation. 

The formal solution to this filtering problem is well known, 
and the frequency response of the optimal Wiener filter is 



1 



(e— ) 



(2) 



The functions and $_ are obtained by the spectral 
factorization of the power spectrum <I> such that 



(3) 



and such that is an analytic function in V without any zero 
inside D and such that is an analytic function outside V 
with <^-{z) ^ for all |z| > 1. The symbol [•]+ in (|3 is 
used to indicate the positive-time part of the sequence whose 
Fourier-transform is inside the brackets. So let / be a given 
transfer function with Fourier coefficients fk, then the positive- 
time part [/]+ is given by f+{z) = Y.kLo fk^^- 

'Throughout this paper, this condition is always implicitly assumed, if a 
function space over the interval [— 7r,7r) is considered 



Solution Q is typical interpreted as a cascade of two filters 



with the two causal filters 



HEie-n = 



(4) 



n + 



(5) 

The whitening filter Hw generates at its output a white noise 
process, whereas the succeeding estimation filter He produces 
the desired minimum-MSB solution x from the whitened 
signal. 

To determine the causal Wiener filter (|2}, we have to solve 
two problems: First, the spectral factorization (|3} of $, and 
secondly the positive time part r+ of the function F := /^^ 
has to be calculated. The formal solutions of both problems 
are well known in the engineering literature (see e.g. [7]). The 
spectral factorization Q is given by 



$+(z) 



and 



=$+(1/^) (6) 



for all z E V and in which F+, given by 



F+{z) ^ exp 



1 

2^ 



(7) 



is called outer function of and ^± are called the spectral 
factors of $. The positive time part of F is by definition 
r+(z) = J2T=Q'ykz''. inserting the Fourier coefficients 7fc 
Q of r gives the closed form solution for all z E T) 



1 pi 
27r ,1^ - 



dr . 



(8) 



It should be noted that r+ exits if T is absolute integrable. 
The spectral factorization has a solution, if $ is absolute inte- 
grable and satisfies the Paley-Wiener condition. Note the very 
similar structure of and (|8}. The solution (|8ji is given by 
the Cauchy transform of F with the kernel X]fe°=o[^^ '^l*^ ~ 
giT^j-giT _ ^-j 'pjjg solution Q is given by the exponential 
function of an integral transformation of the logarithm of 
wherein the kernel is (l + Efclil-^e-*^]*-') /2 = (e" + 
z^jie" — z) and therefore very similar to the kernel of the 
Cauchy transform. Nevertheless, even although both solutions 
have such a similar formal structure, the behavior of both 
solutions is very different, as shown in this paper. 

Therewith, the solution of the Wiener filtering problem is 
formally completely solved. From a practical point of view, 
this solution has to satisfy some desired properties such as 
stability or the possibility to approximate it by an FIR filter 
Subsequently, we study which conditions $ and 5* have 
to satisfy such that the Wiener filter is stable and can be 
approximated by an FIR filter. To this end, the two filters 
Hw and He ait considered separately at first. 

III. Properties of the whitening filter 

In this section, the properties of the pre-whitening filter Hw 
are investigated. Its formal solution is given by 



Hw{z) = [F+{z)\ 



(9) 



in which is defined by Q in terms of the given $. 

First we consider the case that $ satisfies a certain Hoelder 
condition of order a. It turns out that the corresponding Hw 
will be in -4(1?) and satisfies the same Hoelder condition 
as $. So in this case, Hw is a stable filter which can 
be approximated uniformly by an FIR filter Moreover, the 
order a of the Hoelder condition determines how fast the 
approximation error decreases as the degree of approximation 
polynomial increases (cf. Section IVTl . In the second part of 
this section, it is shown that there exist continuous functions 
$ such that Hw is not in ^(P). Therefore it is not possible 
to approximate Hw by an FIR filter in this case. 

A. Smooth spectral densities 

Theorem 1: Let $ e Cq[— 7r,7r) with < a < 1 be a real 
valued function with > S > for all uj E [— tt, tt). Then 
for its outer function F_|_ given by Q holds: e Aai'D) 
and \F+{z)\ > ^ > for all z e V. 

The requirement that $ should be strictly positive is a 
natural condition for our problem, because it ensures that Hw 
always exists. It should be emphasized that for a = 1 the 
above theorem is not valid, in general. 

Before this theorem is proofed, we state the consequences 
of this result concerning the spectral factors $± and the 
whitening filter Hw- Since they are all functions of (see 
(EJ and ©) and since is bounded and strictly positiv in 
V, it can easily be shown that ^± and Hw belongs to same 
space AaiT)) as F+. Thus: 

Corollary 2: Let $ e Cq,[— tt, tt) for < a < 1 be a real 
valued function with > (5 > for all lo E [— tt, tt). Then 
<i>+(z), $_(l/z) and Hw{z) belong to AaiV) and it holds 
|$+(z)| > V5 and \^-{l/z)\ > VS for all \z\ < 1. 

This shows that if <i> is Hoelder continuous, the spectral 
factors <i>± and the whitening filter Hw are stable and can be 
approximated arbitrarily well by FIR filters. Moreover, they 
satisfy the same Hoelder condition as $. 

Proof: (of Theorem [0 From the assumptions on $ 
follows immediately that log<l> e Cq, [—7r,7r). Consider then 
F+ as defined in and write it at z = re*'^ as 



F+{re"^) = exp ( u(re'") + i ■ v{re"^) ) 



(10) 



with 



ZTT 1 — 2r cos((jj — T) + 

The function u is the Poisson integral of log <I> and therefore u 
is harmonic in V and continuous in V with u(e*'^) = log $(w) 
for all u! e [— TT, tt). Now, it can be shown that u satisfies 
in the whole V the same Hoelder condition as log$ and 
therefore u £ Cq,(I?). The function v is the conjugate Poisson 
integral of log<I> or equivalent v is the Poisson integral of 
the Hilbert transform of log<i>. Next, we use a theorem of 
Privalov-Zygmund (see e.g. [5, chapter 3]) which states that 
if a function $ belongs to Ca[— tt, 7r),0 < a < 1, so does 



its Hilbert transform as well. For our case follows that the 
Hilbert transform of log $ belongs to Ca [— tt, tt) and therefore 
with the same arguments as for u: v G Ca{T>). Since u and v 
are harmonic functions in T) and Hilbert transforms of each 
other, the complex function g+ := u + i • u is analytic in T) 
and consequently g+ e Aai'D). The outer function is obtained 
by — exp(g-|_). Of course, F+ is also analytic in V and 
continuous in T) and it is easily verified that F+ satisfies also 
a Hoelder condition of order a such that indeed G Aa{'D). 
Finally, — exp (u) with u given by the Poisson integral 
(tm of log<i>. Because $(tj) > 5 for all uj e [— 7r,7r), it 
follows that log<i>(ijj) > log (5 and since $ is continuous, we 
obtain that |i^+(z)| > (5 for all z e P. ■ 

B. Continuous spectral densities 

Now we show that the assumptions on the spectral density 
made in the last sub-section, can not be weakened, i.e. if 
<!> is continuous but not Hoelder continuous, there exists such 
a $ for which Fj^ belongs not to A{T)). 

Since the given spectral density $ e C+[~tt^ tt) is bounded, 
the outer function Q is analytic and bounded in T). This means 
that e 1-L°°{T>) represents a transfer function of a causal 
and stable filter. Moreover it is well known that every F £ 
H°°{'D) has a radial limit limr-^i F{re^'^) almost everywhere 
(a.e.) in [—tt, tt), and it was shown that the absolute value of 
is continuous in the whole T) with 

lim^^i |F+(re'")| = <P{uj) , Vcj S [-tt,tt) . 

Nevertheless, as the following theorem shows, the outer func- 
tion itself is in general not continuous in V and therefore not 
an element of AiV). 

Theorem 3: To any set E C [—tt, tt) of Lebesgue measure 
zero there exists a function $ G C+ [—tt, tt) such that its 
corresponding outer function F+(re"') has no limit for all 
d e E as r ^ 1. 

In view of the spectral factorization and in respect to the 
properties of the whitening filter (|9} the following Corollary 
is deducted from this theorem. 

Corollary 4: There exist spectral densities $ e C+ [—tt, tt) 
such that there are no functions $+(z) G -4(1?) and 
^-{1/z) G A{'D) which factorize $ according to ^ and 
such that Hw belongs not to A{'D). 

For the proof of Theorem |3] the following lemma, which 
was already proofed in [8], will be necessary. 

Lemma 5: Let E C [—tt, tt) be an arbitrary set of Lebesgue 
measure zero. There exits a function h G C+ [—tt, tt) such that 



lim 



1 



.1 2tt 



h{Lu) 



r sin('i9 — uj) 



1 - 2rcos(i? - w) + H 



dcu = —oo (13) 



for all d G E. 

Proof: (Sketch of proof of Theorem O Let E G [—tt,tt) 
be an arbitrary set of Lebesgue measure zero and let /{u;) := 
exp{h{uj)) for all u; G [—tt,tt), where ft, is a function 
corresponding to the set E as in Lemma |5] Note that the 
so defined / is a continuous and strictly positive function. 



Consider now as defined in (Q and write it at 2: = re^'^ 
as in ilQ\ with u and v as given by (II 1> and ( I12t . Using 
Lemma |5l for v, it can be shown that for all -d ^ E 

limsup 5R{F+(re"')} = $(1?) (14) 
liminf 5R{F+(re"')| = (15) 

And a completely similarly result holds for the imaginary part. 
Since ^ 0, this proofs that both ^{F+} and 

and therefore ^+(19) itself diverges for all -d E E. ■ 
Theorem |3] states that there exists a $ e C+ [— tt, tt) 
such that the corresponding belongs not to ^(P). Since 
the spectral factors $± and the whitening filter Hw are 
determined by (see Q), Corollary 0] follows at once. 
This result has important practical consequences, because even 
if the spectral density $ is continuous, it may not be possible 
to realize the desired filter characteristic of Hw by an FIR 
filter with arbitrarily small error. 

IV. Properties of the estimation filter 

Now the properties of the estimation filter He, given by 
(|5}, are investigated. As we will show, there is a fundamental 
difference between the whitening filter H\y and He- As 
it was shown, Hw is always stable. However, the stability 
of He depends strongly on the smoothness of the given 
spectral densities. If they satisfy no Hoelder condition. He 
may become unstable. First we consider again the case where 
$ and vj/ satisfy a certain Hoelder condition and show that 
He is also a continuous function in V and satisfies the same 
Hoelder condition as the given spectra. In the second part, 
spectra are considered which satisfy no Hoelder condition. It 
turns out that in this case. He is even not bounded, in general. 

Because of the limiting space, the proofs are omitted in 
this section. However, because of the similarity of the two 
operators Q and (|8j, the techniques with that the results can 
be proofed is very similar to that used in Section IIIIl 

A. Smooth spectral densities 

We assume that $ e Ca[-TT,n) with > for 

all u e [— TT, tt). From this assumption follows (see Corol- 
lary |2j that the spectral factor $_(l/z) belongs to ^1^(1?) 
and its absolute value is strictly positive in V. In particular 
$_(e"^) e Ca[-7r,7r) with |$-(e"^)| > Oforallw e [-7r,7r). 
Therefore, also l/$_(e"^) belongs to Cq[— tt, tt). Furthermore 
it is assumed that e Cq[— tt, tt). Since Cq[— tt, tt) is an 
algebra, the quotient F = belongs also to Ca[— tt, tt). 

But to which function space belongs the positive time part of 
r? The answer is given by 

Theorem 6: Let F S Cq[— 7r,7r) with < a < 1. Then the 
positive time part F+ belongs to Aa{T>). 

This theorem shows that as long as both spectra $ and ^ 
satisfy a Hoelder condition. He will be stable and can be 
approximated arbitrarily well by an FIR filter 



B. Continuous spectral densities 

It was shown (Corollary |4j that if $ is only in C[— tt, tt), 
the corresponding spectral factor $_ may not be continuous in 
T). However, we assume that F = 'I'/<I>_ belongs to C[— tt, tt). 
For this case, we have 

Theorem 7: To any set E C [— tt, tt) of Lebesgue mea- 
sure zero there exists a function F G C[— tt, tt) such that 
|F+(re*'')| — > 00 as r ^ 1 for all -d E E and therefore 
||F+||^=oo, i.e. T+in°°{V). 

This theorem shows that if $ and 4* satisfy no Hoelder 
condition, it can not be guaranteed that the estimation filter 
He is stable in the energy sense. 

V. Properties of the Wiener Filter 

Now we consider the properties of the overall Wiener filter 
(0}. They are easily obtained by combining the results of the 
two preceding sections. There, it was shown that if $ and 
are sufficiently smooth, i.e. satisfy a Hoelder condition of order 
a, the whitening filter Hw and the estimation filter He will 
belong to Aa{'D). Under point- wise multiplication, Aa{T>) is 
an algebra. Therefore, also the Wiener filter H = Hw ■ He 
belongs to Aa{T>). This important result is summarized in the 
following corollary: 

Corollary 8: Let £ Ca [— tt, tt) then the transfer func- 
tion H of the corresponding Wiener filter belongs to Aa{T)). 

This corollary says that under the above assumptions on 
$ and 5*, the Wiener filter (|2j is always stable and can be 
approximated arbitrary well by an FIR filter. Moreover, it 
makes a statement on the smoothness of H. As it is shown 
in the next section, from the order a of the Hoelder condition 
the dependency of the approximation error from the FIR 
filter order can be determined. If $ and ^ are not Hoelder 
continuous, the Wiener filter my not be stable. Also this result 
is summarized in a corollary: 

Corollary 9: There exists continuous spectra $ e C[— tt, tt) 
and ^I* e C[— TT, tt) such that the Wiener filter H, given by (|2j, 
is not stable, i.e. such that H ^ n°°{V). 

Of course, because of Corollary |8] at least one of $ and 
satisfies no Hoelder condition in this case. 

Finally, it should be mentioned that all the results, derived 
here for continuous spectral densities which satisfy a Hoelder 
condition of order a (i.e. functions of Ca) can easily be 
extended to spectral densities which are fc-times continuous 
differentiable, and whose fc-th derivate satisfies a Hoelder 
condition of order a (i.e. functions of C^). 

VI. FIR Approximation 

Corollary |8] shows that if $ and ^E* are sufficiently smooth, 
the Wiener filter H belongs to A{'D) and can therefore be 
approximated uniformly by a polynomial. The minimal error 
(measured in the stability norm) which can be achieved by 
approximating H E -4(1?) by an polynomial of degree N 
(Vn denotes the set of all such polynomials) is denoted by 

eN{H)=mipev.\\H-P\\^. 



eAr(iJ) is called the best approximation of H of order N, 
and due to a theorem of Weierstrass, 6^(11) tends to zero as 
N oo for all H e yl(2?). Moreover, for functions which 
satisfy additionally a Hoelder condition of a certain order a, 
the order with that goes to zero is known [5]: 

Lemma 10 (Jackson,1912): Let H g Aai'D) then 
€m{H) < CN^" with a constant C independent on N. And 
if to a given there exists a polynomial approximation of 
degree N such that ejv(i?) < CiV"", then e yla(2?)- 

Thus, as smoother H (i.e. as larger a) as faster decreases 
the approximation error as N is increased. 

The question is now; Which polynomial of degree N 
achieves the best approximation €n{H)1 Or at least, which 
polynomials show the approximation behavior of Lemma [TOf 
To answer this question, we write H e A{'D) as H{z) — 
Sfc^o ^kZ^ and consider the so called Fejer mean of H (see 
e.g. [5, chapter 3]) 



k 

N+l 



Therewith a class of approximation polynomials can be de- 
fined, for which the error decreases as in Lemma [Tol 



Vn{z;H) = 2F2N-i{z;H) - Fn-i{z;H) 



(16) 



They are called the de la Vallee-Poussin means of H{z) and 
the following result is well known. [5, chapter 3] 

Lemma 11: Let Vn be the de la Vallee-Poussin mean of 
H, then \\H-VN\\^<ieN{H). 

For H e AaiTy) follows (together with Lemma [Tot for the 
approximation error that \\H — Vn\\^ < 4CiV^". Moreover 
the approximation error of the de la Vallee-Poussin mean 
V]\[{z, H) is upper bounded by the best approximation e]y{H), 
although it should be noted that Vn{z;H) is a polynomial 
of order 2N — 1. The consequences for our three filters are 
summarized in the following two theorems: 

Theorem 12: Let $, e Ca[-n,T:) and let H be the 
transfer function of the whitening-, of the estimation- (|5}, or of 
the Wiener filter Then, there exists an FIR approximation 
Hn of degree N such that \\H ~ Hn\\^ < CN-" with a 
constant C independent of N. Moreover, this approximation 
velocity is achieved by using de la Vallee-Poussin means 
Vn{z,H) of H as approximation polynomials _ffjv- 

With the previous discussion, this theorem is already 
proofed. For the whitening filter Hw also the reverse is true: 

Theorem 13: Let Hw be the whitening filter (|9j corre- 
sponding to a given spectral density function $ and let Hn 
be an FIR approximation of degree N such that 



\H 



w 



Hn\\^<CN' 



then $ e Cal-TT, n). 

Altogether this shows that the smoothness of $ and is 
directly related to the rate with that the approximation error 
approaches zero. As smoother the densities are, as shorter the 
approximative FIR filter can be. 

However, if the spectral densities are not Hoelder continu- 
ous, the corresponding filters can not be approximated by an 



FIR filter. Moreover, it can be shown by using Theorem|3]and 
Corollary |4] that there exist spectral densities $ € C+ [— 7r,7r) 
such that the best approximation of its outer function F+, using 
an approximation function g E A{'D) from the disk algebra is 



inf 



F. 



\F. 



and that this minimal approximation error is achieved by 
g = 0. Of course this means that can actually not be 
approximated by a function from AiV) and therefore it can 
not be approximated by an FIR filter 

Finally, we mention that the methods of this paper can also 
be used to study the continuous dependency of the Wiener 
filter H on the given data $ and ^E*. First results show that 
every $ G C+[— tt, tt) is a discontinuity point of the spectral 
factorization. This means that for every $ with ||$ — $||oo < e, 
the relative error in the spectral factors ^||$+ — 'I'+lloo may 
becomes arbitrarily large. Similar effects appear also in FIR 
filter design. There, the error grows proportional to log N ■ e 
where N is the FIR filter degree. A detailed investigations will 
be published elsewhere. 

An extension of the investigations of this paper to multiple- 
input multiple-output systems will be given in [9], [10]. 

VII. Conclusions 

For a Wiener filter, it is sufficient that the given spectral 
density functions $ and ^I^ satisfy a Hoelder condition, in order 
that the Wiener filter is always stable and can be uniformly 
approximated by an FIR filter Moreover the smoothness 
of $ and ^I* characterize how fast the approximation error 
decreases as the degree of the approximation filter is increased. 
If $ and vj/ are only continuous, the estimation filter may 
become unstable and it may not be possible to approximate 
the whitening filter by an FIR filter It can be shown that this 
holds also for absolute continuous spectra. 
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